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Developments in design sensitivity analysis (DSA) method have been made using two 
fundamentally different approaches as shown in figure 1 . In the first approach, a 
discretized structural finite element model is used to carry out DSA. There are three 
different methods in the discrete DSA approach: finite difference, semi-analytical, and 
analytical methods. The finite difference method is a popular one due to its simplicity, 
but a serious shortcoming of the method is the uncertainty in the choice of a 
perturbation step size of design variables (ref. 1). In the semi-analytical method, the 
derivatives of stiffness matrix is computed by finite differences (refs. 2-4) whereas in 
the analytical method, the derivatives are obtained analytically. For the shape design 
variable, computation of analytical derivative of stiffness matrix is quite costly (ref. 

1). Because of this, the semi-analytical method is a popular choice in discrete shape DSA 
approach (refs. 3 and 4). However, recently, Barthelemy and Haftka (ref. 5) presented 
that the semi-analytical method can have serious accuracy problems for shape design 
variables in structures modeled by beam, plate, truss, frame, and solid elements. They 
found that accuracy problems occur even for a simple cantilever beam. In the second 
approach, a continuum model of the structure is used to carry out DSA. For shape design 
variable, the material derivative concept of continuum mechanics is used to relate 
variations in structural shape to measures of structural performance (refs. 6-10). 
Using continuum DSA approach, expressions for shape design sensitivity are obtained in 
the form of integrals with integrands written in terms of natural physical quantities 
such as displacements, stresses, strains, and domain shape changes. If exact solutions of 
the continuum equilibrium equations are used to evaluate these continuum design 
sensitivity expressions, the method is called continuum-continuum (C-C) method. On 
the other hand, if the analysis results of the finite element or boundary element methods 
are used to evaluate these terms, the method is called continuum-discrete (C-D) method. 
The analytical method of discrete design sensitivity analysis approach will be called 
discrete-discrete (D-D) method. 


METHODS OF DESIGN SENSITIVITY ANALYSIS 
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The D-D method starts with the finite element matrix equilibrium equation for 
linear structural system as shown in figure 2, where K(b) is the reduced global 
stiffness matrix, z is the reduced displacement vector, F(b) is the external load vector, 
and b is a design variable vector. Differentiating both sides of the matrix equilibrium 
equation with respect to b, a matrix equation for the derivative of displacement vector, 
dz/db, is obtained where the tilde (~) indicates a variable that is to be held constant for 
the process of partial differentiation. If the derivative dz/db is obtained by solving this 
equation, the method is called direct differentiation method. If derivatives of a general 
performance measure are needed, an adjoint variable method can be used (ref. 11). Even 
though the direct differentiation and adjoint variable methods are different in 
computational efficiency depending on situations, they are equivalent in accuracy as Jong 
as consistent computational procedure is used for both methods. For the D-D method, the 
derivative of stiffness matrix is obtained analytically, whereas it is obtained by finite 
differences for the semi-analytical method. The discrete DSA approach is applicable to 
both sizing and shape design variables. For the shape design case, the design variables 
are positions of the finite element grid points. 


DISCRETE DSA APPROACH (DIRECT DIFFERENTIATION METHOD) 

K(b)z = F(b) 

K(b>4r = - -J-(K(b)z) + 
d b 9b 9b 


• Semi-analytical and analytical (D-D) methods 

• For general performance measures, use the adjoint variable method. 

• Accuracy of the direct differentiation and adjoint variable methods are 
equivalent. 

• Discrete approach is applicable to both sizing and shape design variables. 


Figure 2 
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For the continuum approach, using the principle of virtual work, the variational 
equilibrium equation of the structural system can be obtained (ref. 11) as shown in 
figure 3, where a a (-,-) denotes energy bilinear form, i n (-) denotes load linear form, Q 

is the shape of the structure, z is the displacement, z is the kinematically admissible 
virtual displacement, and Z is the space of kinematically admissible virtual 
displacements. Note that an approximate finite element matrix equilibrium equation can 
be obtained by applying the Galerkin method to the variational equilibrium equation for 
an approximate solution. For shape DSA, taking the material derivative of both sides of 
the variational equilibrium equation (refs. 10-12), a variational equation for the 

material derivative z of the displacement is obtained where V is the design velocity field. 

Expressions for a^(z,l)and i^fl) can be obtained for various structural components 

(refs. 10-12). For the C-D method, an approximate finite element matrix equation is 

used to obtain an approximate solution of the second variational equation for z. On the 
other hand, for the C-C method, the analytical solution z of the first variational equation 

is used in the second variational equation to obtain the analytical solution z. As in the D- 
D method, if derivatives of a general performance measure are needed, an adjoint 
variable method can be used (refs. 10-12). The C-C method provides the exact design 
sensitivity of the exact model, whereas the C-D method provides an approximate design 
sensitivity of the exact model. On the other hand, D-D method yields the exact design 
sensitivity of an approximate finite element model, and both the finite difference and 
semi-analytical methods yield approximate design sensitivities of an approximate finite 
element model. 


CONTINUUM SHAPE DSA APPROACH (DIRECT DIFFERENTIATION METHOD) 


a n (z,z) = i n (z), for all zeZ 
a n (z,z) = iy (z) - a^(z,z), for all zeZ 


« FEM equation is an approximate equation of the variational equation. 

* Use material derivative concept of the continuum mechanics for shape DSA. 

* For general performance measures, use the adjoint variable method of DSA. 

* C-C and C-D methods 


Figure 3 
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One question often asked is; "Are the D-D and C-D methods equivalent?" For this 
question, certain conditions have to be given. First, the same discretization (shape 
function) used for the finite element analysis method must be used to evaluate the 
continuum design sensitivity results. Second, exact integrations (instead of numerical 
integrations) must be carried out for all integrations used for generation of stiffness 
matrix and evaluation of continuum design sensitivity expressions. The third condition to be 
met is that the exact solutions (not a numerical solution) of the finite element matrix 
equation and adjoint equation are used to compare two methods. The fourth condition is that 
movement of the finite element grid points for shape design change in the D-D method 
must be consistent with the parameterization method used for the design velocity field of 
the C-D method. For the sizing design variable, it is shown in reference 1 1 that the D-D 
and C-D methods are equivalent under the conditions given in figure 4 using a beam 
structural component. It has also been argued that the D-D and C-D methods are 
equivalent for shape design variable under the conditions given in figure 4 (refs. 13 and 
14). One point to note is that these four conditions are not easy to satisfy; in many cases, 
numerical integrations are used and exact solutions of the finite element matrix 
equations cannot be obtained. In this paper, equivalence study of D-D and C-D methods 
for shape design variables is carried out under the conditions given in figure 4. To carry 
out equivalence study of the D-D and C-D method, two simple structural components, a 
truss and a cantilever beam, are used. The shape DSA results of the D-D and C-D methods 
derived in the published literature are cited and used here without being derived in this 
paper. 


ARE THE D-D AND C-D METHODS EQUIVALENT? 


Equivalence study under the following conditions: 

1 The same shape function used for FEA must be used to evaluate the continuum 
DSA results. 

2 Exact integrations must be used to generate the stiffness matrix and evaluate 
the continuum DSA results. 

3 Exact solutions of the finite element and adjoint matrix equations are used to 
compare two DSA methods. 

4 Movement of FE grid points for the D-D method must be consistent with the 
parameterization of the design velocity field for the C-D method. 


Figure 4 
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In figure 5, the results of equivalence study of the D-D and C-D methods for shape 
design sensitivity are presented using a simple truss with one end fixed. The truss has a 
uniform cross-sectional area A and its length is i which is the shape design variable. 
Three loading cases, a point load p at the tip, a uniformly distributed load f, and a 
linearly varying load qx/i, are considered as shown in figure 5. For each loading case, 
linear and quadratic shape functions are used for finite element models. For the linear 
shape function, two element model is used whereas for the quadratic shape function, one 
element model is used. For the equivalence study, design sensitivities of the nodal 
displacements are considered using the adjoint variable method. In figure 5, 'same' 
denotes that the D-D and C-D methods yield the same result and 'not' denotes that the two 
methods do not yield the same result. Details of the equivalence study results are given in 
the following figures. 


RESULTS OF EQUIVALENCE STUDY OF D-D AND C-D METHODS FOR TRUSS 
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Note that for all three parameterizatlons of design velocity, the perturbation of length of 
the truss is 64 at the tip. Moreover, for the quadratic and cubic design velocities, once 
C], 1-1 ,2,3, are fixed, then the only one shape design variable is the length X. The 
movement of the finite element grid points for shape design changes in the D-D method 
must be consistent with these parameterization methods. For the D-D method, the shape 
design variables are the positions and b 2 of the nodal points. If the present design Is 

bt-/l/2 and b 2 -i, then V(i/2)a6bi -6,1/2 and V(/6)a6b 2 -64 for the linear velocity, 
V(i/2)a6bi-Ci and V(i)a6b 2 -6i for the quadratic velocity, and V(i/2)a6bi- 
(9 i 2 + 9 e 3 * 6i)/1 6 and V(i)a6b 2 «6i for the cubic velocity. 


PARAMETERIZATIONS OF THE DESIGN VELOCITY V(x) 



0 Jt / 2 X 0 XI 2 X o 4/32^/3 4 

Linear Quadratic Cubic 


V^x) = x6i/i V 2 (x) - a 2 x 2 + o^x V 3 (x) - p 3 x 3 + p 2 x 2 + pjx 

Figure 6 
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The first case of equivalence study is the truss with the point load p at the tip. For 
this, the finite element matrix equation, using linear shape function, is given in figure 7 
where the stiffness matrix depends on the shape design variables bj, i=1,2. The finite 

element matrix equation gives the solutions z-\ sz(i/2)=p/£/2EA and Z 2 *z(i)=pi/EA at 
the present design b-|=i/2 and b 2 =i. Thus z(x)=px/EA which is the exact solution of 
the truss with the point load p at the tip. If the design sensitivities of displacements at 
two nodal points, z-\ and z 2 , are desired, the adjoint equations are given in figure 7 with 

the adjoint solutions XJ(x)=x/EA for 0<x<X/2, X 1 (x)=/8/2EA for X/2<x<X, and 
X 2 (x)=x/EA, respectively. These adjoint solutions are also exact. 


FIRST CASE: TRUSS WITH THE POINT LOAD p AT THE TIP 
(Linear Shape Function) 
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Using the D-D method, design sensitivities for z-| and Z2 are z' -p6i/2EA and 
z^=p6i/EA, respectively, for the linear velocity as shown in figure 8. On the other 
hand, if the quadratic velocity is used, then z^=p£i/EA and z^=p6i/EA. Also for the 
cubic velocity, the D-D method yields z^=p(9E 2 +9£3-6 i)/16EA and z^=p6i/EA. Now, 

using the C-D method, the design sensitivity expression is obtained as 

x 



Using the finite element analyses results and the linear velocity in this design 
sensitivity expression, the C-D method gives z^=p6i/2EA and z^=p6i/EA which are 

the same as the results of the D-D method. Moreover, the design sensitivity expression 
yields z^=p£i/EA and z^=p6i/EA for the quadratic velocity and z' i =p(9£2+9€3- 

6i)/16EA and z^=p6i/EA for the cubic velocity which are the same as the results of 

the D-D method. Thus, when the linear shape function is used for finite element model of 
the truss with the point load p, the D-D and C-D methods are equivalent for all 
parameterizations of velocity considered as indicated in the second column of figure 5. 
One point to emphasize in this case is that the original and adjoint responses of finite 
element models are the exact solutions of the truss with the point load. Note that the 
design sensitivity z^=p6i/EA is independent of the parameterizations of velocity for the 

C-D method. 


DESIGN SENSITIVITY OF NODAL DISPLACEMENTS (First Case) 


* D-D and C-D methods yield the same result for all parameterizations of 
velocity. 


Linear Velocity 
Quadratic Velocity 
Cubic Velocity 


z^=p6i/2EA 

z’=p£i/EA 

z’=p(9£ 2 +9E 3 -6i)/16EA 


Z 2 =p6i/EA 

Z 2 =p6i/EA 

Z 2 =p6i/EA 


Figure 8 



The second case of study is the truss with the uniformly distributed load f along the 
truss. For this, using the quadratic shape function, the finite element matrix equation is 
obtained as given in figure 9. The solutions of the finite element matrix equation are 
zi»z(i/2)=3fi 2 /8EA and Z2=z(i)=fi 2 /2EA at the present design bfi/2 and b 2 -i. 

Thus z(x)=fx(-x+2i )/2EA which is the exact solution of the truss with the uniformly 
distributed load f. if the design sensitivities of zt and Z 2 are desired, the adjoint 

equations are given in figure 9 with the adjoint solutions X^(x)=(-3x z /4X+5x/4)/EA 
and X 2 (x)=x/EA, respectively. The adjoint solution X 2 (x) is the same as in the linear 
shape function case which is the exact solution. On the other hand, the adjoint solution 
X 1 (x) is different from the linear shape function case and not exact. 


SECOND CASE: TRUSS WITH UNIFORMLY DISTRIBUTED LOAD f 
(Quadratic Shape Function) 
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Using the D-D method, design sensitivities for z-\ and z 2 are z^=3fi6i/4EA and 
z^=fi6i/EA, respectively, for the linear velocity as shown in figure 10. On the other 
hand, if the quadratic velocity is used, then z =fi(6i+e 1 )/2EA and z'=fi6i/EA. Also, 
for the cubic velocity, the D-D method yields z'^-fi (156i+9€ 2 +9£3)/32EA. Now, 

using the C-D method, the design sensitivity expression is obtained as 

x 

Zj' = f ( f X. 1 + EA z x xj< ) V x dx, i = 1 ,2 

Jo 

Using the finite element analyses results and the linear velocity in this expression, the 

C-D method gives z =3fi6i/4EA and z =fi6i/EA which are the same as the results of 
1 2 

the D-D method. Also, using the finite element analyses results and the quadratic 
velocity in the design sensitivity expression, the C-D method gives z^=fi(6/8 + £i)/2EA 

and z =fi6i/EA. These are the same as the results of the D-D method. However, the 
2 

design sensitivity expression yields zj*fi(426i+36£ 2 +9£ 3 )/80EA for the cubic 

velocity which is different from the result of the D-D method. Hence, it can be concluded 
that the D-D and C-D methods are not equivalent in the second case of study. Notice that 
the sensitivity results of the D-D method are the same as those of the C-D method up to 
the linear velocity when the linear shape function is used and up to the quadratic velocity 
when the quadratic shape function is used. Thus, the second case indicates the D-D and C- 
D methods might be equivalent under an additional condition that the shape function used 
in the finite element model is isoparametric with the discretization polynomial of the 
design velocity. However, this is not true as the results of the next case of study indicate. 


DESIGN SENSITIVITY OF NODAL DISPLACEMENTS (Second Case) 


D-D 


C-D 


Linear Velocity 


z =3f/£ 6/6/4EA 
1 

z =fi6i/EA 
2 


z'-3fJ6J/4EA 

z =fi64/EA 
2 


Quadratic Velocity 


zJ=fi(6i+£ 1 )/2EA 

z’ =fi6i/EA 
2 


z^U(6i+£ 1 )/2EA 

z =fi6i/EA 
2 


Cubic Velocity 


z^=fi(156i+9£ 2 +9£ 3 )/32EA 


z’=fi(426i+36£ 2 +9£ 3 )/80EA 


Figure 10 



The third case of study is the truss with the linearly varying load qx/i along the 
truss. Before carrying out design sensitivity computation, dependency of the external 
load on the shape design has to be defined as shown in figure 1 1 . That is, as the length of 
the truss changes, the external load will maintain the form of qx/i. For this, using the 
quadratic shape function, the finite element matrix equation is given in figure 1 1 . The 
matrix equation gives the solutions Zi=z(i/2)=11qi 2 /48EA and z 2 =z(i)=qi 2 /3EA at 
the present design b^i/2 and b 2 =i. Thus z(x)=qx(-3x+7i)/12EA which is not the 
exact solution of the truss with the linearly varying load. The same adjoint equations that 
are given in figure 9 are applicable in this case with the solutions (x)=(-3x 2 /4A 
+5x/4)/EA and A, 2 (x)=x/EA, respectively. As mentioned before, the adjoint solution 
A, 1 (x) is not exact, whereas A 2 (x) is exact. 


THIRD CASE: TRUSS WITH LINEARLY VARYING LOAD qx/4 
(Quadratic Shape Function) 
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Using the D-D method, design sensitivities for z-\ and Z2 are z^=11qi6i/16EA and 

z^=qi6i/EA, respectively, for the linear velocity as shown in figure 12. On the other 

hand, if the quadratic velocity is used, then z^=qi(256/8+16£i)/48EA. Also for the 

cubic velocity, the D-D method yields z^=qi6/8/EA. Now, using the C-D method, the 

design sensitivity expression is obtained as 
x 




x+EAz^V, 


dx, i = 1 ,2 


Using the finite element analyses results and the linear velocity in this expression, the 
C-D method gives z^=1 1qi6i/16EA and z^=qi6i/EA which are the same as the results 

of the D-D method. However, using the finite element analyses results and the quadratic 
velocity in the design sensitivity expression, the the C-D method gives 


z’=qi(196i+292£ 1 )/240EA, whereas it yields z^=qi(2496i+27€ 2 -27€ 3 )/240EA 

for the cubic velocity. These are not the same as the results of the D-D method. Thus the 
D-D and C-D methods are not equivalent for the truss with a linearly varying load. Based 
on the equivalence study of truss problem, the D-D and C-D methods are possibly 
equivalent only for linear velocity. If this is the case, then both methods will give the 
exact design sensitivity information of the finite element analysis results that may not 
be acceptable at all. This is the situation for the fillet problem in reference 15 that the 
design sensitivity results of the C-D method agrees up to 5 to 6 digits with the finite 
difference even though the finite element model using constant stress triangular element 
does not provide accurate analysis result. On the other hand, when automatic regridding 
methods are employed for shape optimal design (refs. 16 and 17), parameterizations of 
the design velocity field cannot be limited to be only linear functions. 


DESIGN SENSITIVITY OF NODAL DISPLACEMENTS (Third Case) 
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The results of analytical equivalence study for a simple cantilever beam with 
moment of inertia I and length Jt are given in figure 13. Like the truss problem, three 
lateral loading cases shown in figure 13 are considered. For all loading cases, Hermite 
cubic shape functions are used for the finite element model with one element. Also, for 
the design velocity V(x), the same linear and quadratic parameterizations as in the truss 
problem are used. In addition to these, Hermitian parameterization of the velocity is 
used. That is, if the beam is fixed at x=0 and changes its length by 6 Z at x=i and the 
slope of the velocity is zero at x=0 and 0 at x=i, then the parameterization of the 

velocity is V 4 (x) as shown in figure 13 where y and y are given in terms of i, 6i, and 

3 2 

6. For the equivalence study, design sensitivity of the tip displacement is considered. The 
results of equivalence study are summarized in figure 13. As in the truss case, the finite 
element model for the beam with the point load p at the tip yields the exact solutions of 
the original and adjoint structures. Hence the D-D and C-D methods give the same design 
sensitivity results for all parameterizations of velocity as shown in figure 13. 


RESULTS OF EQUIVALENCE STUDY OF D-D AND C-D METHODS FOR BEAM 
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For the beam with the linearly varying load qx/i along the beam, the finite element 
matrix equation is given in figure 14 with the solutions z t =1 1qi 4 /120EI and 
z 2 =qi 3 /8EI at the present design b=i. Thus z(x)=qi(18/6x 2 -7x 3 )/120EI which is an 
approximate solution of the beam with linearly varying load. For the design sensitivity 
of z-), the solution of the adjoint equation given in figure 14 can be used. Using the D-D 

method, design sensitivity for z-| is z^ = 1 1qi 3 6i/24EI for all parameterizations of 

velocity. For the C-D method, the design sensitivity expression is 

x 

z i = | ^ t 3 z xx A. xx V x + ( z x A. xx + z xx A, x ) V xx ] + ^-3. jx.V + j?iV x j dx 

Using the finite element analyses results in this design sensitivity expression, the C-D 
method yields z^=1 Iqi 3 6i/24EI for the linear velocity which is the same as the result 

of the D-D method. However, the design sensitivity expression yields z’ i =qi 3 (1 71 6 i- 

1 2€ i )/360 El for the quadratic velocity and z^=qi 3 (1 1 946 i-9i 0)/252OEI for the 

Hermitian velocity which are not the same as the results of the D-D method. Thus the D- 
D method and C-D methods are not equivalent for the beam with linearly distributed load 
as indicated in figure 13. Based on the equivalence study of the beam problem, the D-D 
and C-D methods are possibly equivalent only for linear velocity. 

DESIGN SENSITIVITY OF NODAL DISPLACEMENT FOR BEAM 
(Linearly Varying Load q x/Jt) 
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Next, a numerical study is carried out for the C-D method using the cantilever beam 
with the uniformly distributed load to see effect of accuracy of the finite element 
analysis results on accuracy of the design sensitivity informations obtained. The finite 
element models with 1 , 2, and 20 elements are considered for numerical study. Node 
numbering for all finite element models starts at the clamped end of the beam and the 
node number of free end of the beam is (m+1) where m is the number of elements in the 
model. The beam is 60 in. long and has a uniform rectangular cross-section of 0.5 in. 
high and 0.25 in. wide. Young's modulus, Poisson's ratio, and uniformly distributed load 

are E=30x10 6 psi, v=0.3, and f=0.5 ib/in., respectively. Finite element analysis is 
carried out using ANSYS finite element STIF4. Three parameterizations of velocity with 
1% perturbation of the length i=60 in. of the beam are used for numerical study as 
shown in figure 15. Once the solutions of the original and adjoint structural system are 
obtained using ANSYS, the continuum design sensitivity expression is numerically 
integrated using three points Gauss quadrature. 


CONTINUUM-DISCRETE METHOD FOR A CANTILEVER BEAM 
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Z =60 in., E=30x 1 0 6 psi, v =0.3, f=0.5 ib/in. 

ANSYS STIF4 



Parameterizations of Velocity for Numerical Study of the C-D Method 


Case 

Design 

Velocity Type 

Parameter Values 

A 

Linear 

=0.6 in. 

B 

Quadratic 

6/6 =0.6 in. and z^= 10 in. 

C 

Hermitian 

6 Z =0.6 in. and 0 = -0.3 


Figure 15 
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To check accuracy of the design sensitivity obtained, the results are compared with 
the results obtained by finite difference as shown in figures 16 and 17. In these figures, 
z(i-6i) and z(i+6i) are the displacements of selected nodal points for perturbed 
backward and forward designs, respectively, Az*z(i+6i)-z(i-6i) is the finite 
difference, and z' is the difference predicted by the design sensitivity. The ratio of z' and 
Az times 1 00 can be used as a measure of accuracy of the design sensitivity. In figures 
16 and 17, for all finite element models, the case A with linear velocity yields excellent 
agreement between the design sensitivity z' and the finite difference Az. This confirms 
with the results of analytic study that the D-D and C-D methods may be equivalent for 
linear velocity. On the other hand, for one element model, the design sensitivity z' and 

the finite difference Az do not agree at all for other parameterizations (cases B and C) of 
velocity as can be seen in figure 16. For cases B and C, the agreements improve 
substantially for two elements model. 


COMPARISON OF DESIGN SENSITIVITY OF THE C-D METHOD 
One Element Model 



Node 




(Z’/AZ 

Case 

No. z(i-5,$) 

z(i+8i) 

AZ 

z' 

xl 00)% 

A 

2 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.41 471 E-01 

100.0 

B 

2 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

-0.47926E-01 

-115.6 

C 

2 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.64970E-01 

156.6 


Two Elements Model 
Node 

Case No. z(i-5 X) 

z(i+6i) 

AZ 

z' 

(Z'/AZ 

x100)% 

A 

2 

0.35273E+00 

0.3821 0E+00 

0.1 4689E-01 

0.1 4688E-01 

100.0 


3 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.41471 E-01 

100.0 

B 

2 

0.1 7939E+00 

0.59469E+00 

0.20765E+00 

0.20744E+00 

99~9 


3 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.35880E-01 

86.5 

C 

~2~ 

0.30947E+00 

0.42955E+00 

0.60040E-01 

0.611 26E-01 

101.8 


3 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.42939E-01 

103.5 


Figure 16 



On the other hand, for twenty elements model, agreements become excellent as 
shown in figure 17. This confirms the fact that accurate design sensitivity informations 
can be obtained as long as accurate finite element analysis results are used for the C-D 
methods. This fact is not the case for the semi-analytic method, as demonstrated by 
Barthelemy and Haftka (ref. 5). They found that the design sensitivity error of the 
semi-analytic method is proportional to the square of the number of elements. This is 
completely opposite behavior from the C-D method since the design sensitivity error 
increases very rapidly as the finite element analysis results of the original structure 
become more accurate. As demonstrated in figures 1 6 and 1 7, an essential advantage that 
may accrue in the C-D method is associated with the ability to identify the effect of 
numerical error associated with finite element analysis results. That is, if disagreement 
arises between the design sensitivity of the C-D method and the finite difference, then 
error has crept into the finite element approximation. If the D-D method is used, in 
which the structure is discretized and the design variables are imbedded into the 
stiffness matrix, then any error inherent in the finite element model is consistently 
parameterized and will never be reported to the user. Therefore, precise design 
sensitivity coefficients of the matrix model of the structure are obtained without 
realizing that there may be substantial inherent error in the original model. On the 
other hand, the C-D method can be used to obtain a warning that approximation error is 
creeping into the finite element model. 


COMPARISON OF DESIGN SENSITIVITY OF THE C-D METHOD (Cont) 
Twenty Elements Model 

Node (z'/az 


Case No. z(i-5i) 

z(Jt+SJt) 

A Z 

z' 

x100)% 

A 2 

0.481 58E-02 

0.521 69E-02 

0.20055E-03 

0.20053E-03 

100.0 

6 

0.1 0504E+00 

0.1 1 379E+00 

0.43744E-02 

0.43739E-02 

100.0 

1 1 

0.35273E+00 

0.3821 0E+00 

0.1 4689E-01 

0.1 4688E-01 

100.0 

1 6 

0.66525E+00 

0.72066E+00 

0.27704E-01 

0.27701 E-01 

100.0 

21 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.41471 E-01 

100.0 

B 2 

0.70800E-03 

0.1 3220E-01 

0.62558E-02 

0.62566E-02 

100.0 

6 

0.29727E-01 

0.22933E+00 

0.99800E-01 

0.1 01 27E+00 

101 .5 

1 1 

0.1 7939E+00 

0.59469E+00 

0.20765E+00 

0.21 023E+00 

101 .2 

1 6 

0.50771 E+00 

0.89205E+00 

0.1 921 7E+00 

0.1 9269E+00 

100.3 

21 

0.99593E+00 

0.1 0789E+01 

0.41 476E-01 

0.41 467E-01 

100.0 

C 2 

0.4761 4E-02 

0.52749E-02 

0.25676E-03 

0.25673E-03 

100.0 

6 

0.95061 E-01 

0.1 2480E+00 

0.1 4869E-01 

0.1 4863E-01 

100.0 

1 1 

0.30947E+00 

0.42955E+00 

0.60040E-01 

0.60046E-01 

100.0 

1 6 

0.60858E+00 

0.781 34E+00 

0.86383E-01 

0.86384E-01 

100.0 

21 

0.99593E+00 

0.10789E+01 

0.41 476E-01 

0.41 470E-01 

100.0 


Figure 17 
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